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NONLINEAR  ELECTROMAGNETIC  SCATTERING  FROM  A  DIODE  CIRCUIT 


I.  INTRODUCTION 


The  problem  we  treat  is  that  of  the  electromagnetic 
scattering  by  a  diode  circuit.  This  is  an  inherently 
nonlinear  problem:  indeed,  it  is  the  nonlinear  effects  that 
are  of  interest.  Previous  authors  ( Franceschetti  and  Pinto, 
1980,  Landt  et  al.,  1983,  Liu  and  Tesche,  1976,  Kanda,  1980, 
Sarkar  and  Weiner,  1976)  have  focused  on  scattering  by  a 
nonl inear ly  loaded  antenna.  This  effort  extends  this  previous 
work  by  considering  the  scattering  from  an  active  circuit. 
The  focus  here  is  on  the  interaction  between  the  incident 
fields  and  the  fields  due  to  the  sources  in  the  circuit.  This 
interaction  results  in  a  mixing  of  the  frequency  components  of 
the  incident  fields  and  the  internal  source,  an  effect 
characteristic  of  the  nonlinearity. 

Our  approach  uses  the  Volterra  series  method  for 
approximating  the  solution  of  a  nonlinear  differential 
equation  (Sarkar  and  Weiner,  1976,  Graham  and  Ehrman,  1973, 
Schetzen,  1980).  The  series  solution  is  truncated  at  the 
third  order.  It  is  thus  a  perturbation  method,  where  the 
nonlinear  effects  are  regarded  as  corrections  to  the  more 
dominant  linearities.  Our  solution  is  strictly  applicable  to 
the  mildly  nonlinear  case,  but  it  also  sheds  light  on  the 
behavior  of  the  more  general  case. 


II.  FORMULATION 


We  will  now  formulate  the  circuit  model  and  derive  the 
differential  equation.  The  circuit  we  treat  is  shown  in 
Figure  1.  The  circuit  has  a  sinusoidal  voltage  source  v(t)/  a 
pure  load  resistance  RLand  a  diode.  The  antenna  is  a  simple 
dipole  and  generates  an  additional  voltage  in  the  circuit  due 
to  reception  of  some  incident  RF  signal.  We  will  focus  mainly 
on  the  calculation  of  the  total  scattered  power  at  various 
frequencies  and  the  voltage  across  the  diode. 

Following  Kanda  (1980),  we  will  approximate  the  diode  by 
a  linear  capacitor  C  and  a  nonlinear  resistance  R  in  parallel. 
The  current  through  this  resistor  obeys  the  characteristic 
diode  relation: 

yt)  =  i^e“^i)  (1) 

where  i^  is  the  saturation  current  and  is  a  constant  depending 
on  temperature  and  is  the  voltage  across  the  diode. 
Following  Liu  and  Tesche  (1976),  we  will  represent  the  antenna 
by  its  Norton  equivalent  circuit,  i.e.,  by  a  current  source, 
igj^t)  and  an  admittance  Yin  in  parallel  (where  isc(t)  is  the 
short  circuit  antenna  current  and  Yin  the  antenna  input 
admittance).  We  thus  arrive  at  the  equivalent  circuit 
depicted  in  Figure  2. 

Using  Kirchoff's  current  law  at  the  3  nodes  results  in 
the  differential  equation 

i^t)  =  isc(t)  +  v(t)/RL-[l/R,+Y(D)+CD]  Vjj(t)  (2) 
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where  IrIs  the  current  through  the  nonlinear  resistor,  D  is 
the  differential  operator  d/dt  and  Y(D)  is  a  linear 
differential  operator  due  to  the  antenna  admittance. 
Substituting  from  (1)  produces 

1)  =  isd-  v/Rj_-  [l/Ri^CD+Y(D)]Vjj  (3) 

This  is  the  equation  we  seek  to  solve:  the  solution  will 
yield  '\^(t)  in  terms  of  the  driving  term  ( ig^i  t)+v{  t)/RL  ) .  The 
nonlinearity  of  the  equation  is  due  to  the  exponent.  A 
similar,  but  simpler  equation  is  treated  by  Kanda  (ibid).  The 
nonlinearity  can  be  transformed  into  a  form  which  is  more 
useful  for  our  method  of  solution  by  using  the  series 
expansion  of  the  exponent: 

(4) 

n=0 

Substitution  of  this  into  the  differential  equation  yields 
after  rearrangemment  of  terms: 

(isc+v/Rj)  =  (isa  +  l/R^^^+  C-^+  Y(D)Vj+  igSt^  (5) 

+  •  •  • 

In  this  form,  the  driving  term  is  apparent,  as  well  as  the 
linear  and  nonlinear  parts  of  the  equation. 
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III.  VOLTERRA  SERIES  SOLUTION; 


We  will  now  use  the  Volterra  series  method  (Sarkar  and 
Weiner/  1976,  Graham  and  Ehrman,  1973,  Schetzen,  1980)  to 
derive  an  approximate  solution  to  (5).  The  method  will  utilize 
a  series  of  third  order  and  requires  the  approximation  of  the 
diode  characteristic  equation: 

=„v (6) 

Clearly  this  approximation  of  i  by  a  cubic  is  best  in  the 
vicinity  of  v  =  0,  i.e.,  when 

|av|«:i  (7) 

At  room  temperature,  is  of  the  order  of  38  (volts)  ,  so  our 
solution  will  be  strictly  applicable  for  the  domain 

jv|<::0.026  volts  (8) 

To  simplify  the  following,  let  us  write  the  approximation 
to  (5)  as; 

Y(D)[v(t)]+  CD[v(t)]+  (a, +l/RL)v(t)  +aj^t)  +  (9) 

=  g(t) 

where  we  have  dropped  the  subscript  on  v(t),  the  driving  term 
is  written  as  g(t )  (=isc  +  v/R,)  and  the  diode  expansion 
coefficients  are  a^ =  aj  =  i^2  and  33=  The  Volterra 

series  expansion  for  v(t)  to  third  order  is; 


v(t)  =  v,(t)  +  \^(t)  +  vJt) 


(10) 


vhere 

Vi(  t)  =  j  h/T,)g(t-7i,  dT|  (11) 

’^(t)  =  l^ii/r2)g(t-i^)g(t-'^)dijd'^ 

V3(t)  =  iJiliiTi,'5/T5)g(t~r^)g(t-i5)g(t-i5)dTjd'^r3 

where  all  integrals  are  over  the  range  (-oqoo)  and  where 
hp  ( T^/ . .  .ly  is  the  p***  order  impulse  response  due  to  the  p***  order 
nonlinearity.  Thus/  v^Ct)  is  the  linear  part  of  the  response, 
■\^(t)  is  the  quadratic  and  v^t)  is  the  cubic.  Instead  of 
impulse  responses,  we  may  write  Vp(t)  in  terms  of  transfer 
functions : 

v^(t)  =  l/27r  f  H^(w)G(w)e^w  (12) 

V2(t)  =  l/(27r)^j|l^wi,v^G(wpG(v^e'*^%'WjdW2 

V3(t)  =  l/(27rf  jl^y;wj,>^v^G(w,)G(>^G(>^e^**^fT^*d'h^vgiw3 

Now,  let  g(t)  be  a  sum  of  sinusoids: 

g(t)  =  Ssije''^'*  (13) 

i=i 

Then 

G(w)  =  27r23g|d(w-Wi)  (14) 

is1 

It  can  be  shown  that  the  transfer  functions  are  invariant 
under  a  permutation  of  frequencies  (Schetzen,  1980).  Using 

this  property,  substitution  of  (14)  into  (12)  yields 

Vi(t)  =  EH^(w,)e<'^*  (15) 

V2(  t )  =  z  E  i:  g,  g,  Hjf w„  w,)e*'^‘ 

•  A 

V3(  t )  =  E  g,3H^  V„  w„  w,)e^=»'*i‘  +  3  E  E  91*91  Wj,  Wi,  Wj)e'‘*'*'’"'*'‘“  + 

1 

E  L  p,  9)  9k  m  w,/  w,, 


where  all  sums  are  from  1  to  n.  The  solution  of  our  differ¬ 
ential  equation  (9)  has  thus  been  reduced  to  determining  the 
transfer  functions.  This  is  straightforward  but  algebraically 
tedious  and  will  only  be  outlined  here  (see,  e.g.,  Sarkar  and 
Weiner,  1976,  Graham  and  Ehrman,  1973,  Franceschetti  and 
Pinto,  1980).  Using  (10)  in  (9)  produces 


I^Y  ( D )  +CD+a,+  1/I\j  ( =  g  ( t )  - a^  ^ 


(16) 


With  the  assumption  of  (13),  i.e.,  that  the  driving  term  is  a 
sum  of  sinusoids,  (16)  can  be  solved.  First  we  expand  the 
square  and  cube  in  (16)  and  then  substitute  for  v|,  v^  and  v^  , 
using  (15)  and  for  g(t)  using  (13).  Now  the  effect  of  the 
differential  operators  on  6*'**  is  carried  out.  Terms  are 
grouped  by  frequency  and  order  and  equated  to  zero.  Thus, 
equating  terms  in  yields  the  first  order  transfer 
function; 


“  Y(w)  + jwC+a(*-l/lRL 

Grouping  the  terms  in  6  '  '  produces 


(17) 


W,,  Wj)  =  -  ^2  Hj(  w,)  Wj)  Wj4-Wj) 


(18) 


Grouping  the  terms  in  and  using  (17)  and  (18)  produces 


r^Wj,Wj,w^  =  H(w,+Wj+v^pH(Wi)H(Wj)H(vij^  X 

^2/3a2^  ^H(wj+vj^+H(wj+\iij^H(H+>^)j  -aoj 


(19) 


This  outline  of  the  derivation  of  (18)  and  (19)  has  implicitly 
used  the  property  that  the  p*”  order  transfer  function  depends 
only  on  transfer  functions  of  lower  order.  Equations  (17), 
(18)  and  (19),  together  with  (15),  represent  the  formal 
solution  to  our  equation.  The  presence  of  harmonics  as  well 
as  mixed  frequency  components  is  demonstrated  by  the 
expressions  in  (15).  The  magnitude  of  these  components  can  be 
calculated  from  (17),  (18)  and  (19). 


IV.  SIMULATION; 


We  will  now  use  the  derived  solution  to  simulate  the 
behavior  of  the  circuit  depicted  in  Figure  1.  To  do  so,  we 
must  specify  the  antenna  characterics,  more  specifically,  Y(w) 
and  is^/t).  We  will  restrict  ourselves  to  the  case  where 
plane  waves  are  incident  broadside  to  the  antenna.  This 
simplifies  the  calculation  of  received  current  considerably 
without  obscuring  the  physics  of  the  Volterra  series  solution. 
For  broadside  incidence  of  a  plane  wave  of  frequency  w  and 
amplitude  Eo  on  a  dipole  antenna,  we  have  (Jordan  and  Balmain, 
1968)  ; 


isc(v) 


(20) 


where  the  effective  length.  Left  ,  is  given  by 


T  -  2(1  -  coa  K^) 
"  K  sin  K£ 


(21) 


where  K  =  2Tr/\  =  w/c.  The  calculation  of  the  antenna 
impedance  (or  its  reciprocal,  the  admittance)  is  complicated 
by  the  frequency  range  of  interest;  we  desire  frequencies 
well  below  ,  as  well  as  above,  resonance.  The  method  adopted 
uses  two  different  methods  of  calculaltion.  For  wavelengths 
greater  than  20  times  the  antenna  half-length  (  ,  we  use  the 
short  dipole  approximation  (Tai,  1961); 


Z(w)  =  20K^^-  j  120/Kf  [fn(2f/a)-ll 


(22) 


where  (a)  is  the  antenna  radius.  For  wavelengths  shorter  than 
20£,  we  use  Schelkunof f ' s  method  (Jordan  and  Balmain,  ibid); 


where  Zo=  120  [in( 2//a)-lj 

M  =  6o[^n(2kO-Ci(2k/)+8-l+cos(2k/)] 

N  =  6ofsi(2k/)-sin(2k/)] 

R  =  60r»+in(2k/)-Ci(2k/)] 

+  30|^»+in(ki)-2Ci(2kj2)+Ci(4k^)|  cos(2kj2) 

+  30  [si(4k.C)-2Si(2ki)jsin(2ki) 

X  =  60Si(2kJ)  +  30(Ci(4kj2)  -jCnCkjJ)  -»|sin(2kj2) 


-  30Si  (4kfi  )cos(2kj2) 


and  where  V  is  Euler's  constant  and  Si(x)  and  Ci(x)  are  the 
sine  and  cosine  integrals  respectively: 


»  =  0.5772157 

si(x) 

Ci(x) 


The  use  of  Schelkunoff s  expressions  is  motivated  by  the 
desire  to  have  the  impedance  not  only  at  the  half-wave  and 
full-wave  resonances,  but  also  at  higher  resonances.  In  this 
simulation,  we  will  use  the  following  values: 


H  =  16cm 


a  =  1 . 6mm 


In  Figure  3,  the  antenna  impedance,  using  these  values,  is 
shown  for  frequency  range  1  KHz  to  1  GHz.  The  numerical  data 
from  which  the  figures  were  constructed  showed  that  the 
half-wave  resonance  occurs  at  approximately  450  MHz  and  the 
full-wave  resonance  at  about  820  MHz,  i.e.,  when  Im(Z)  =  0. 


The  calculation  of  the  total  scattered  power  at  frequency 
w  is  easily  done  if  the  voltage  across  the  diode  is  known. 
From  Figure  2,  it  is  apparent  the  voltage  across  the  antenna 
admittance  is  the  same  as  across  the  diode.  Hence,  we  have 
for  total  scattered  power 


P(w)  =  |iiw)|  Ra  (w) 

=  1  V(/w)Y(w)*Ra{w) 
=  |'v^/w)/Z(w)  *Ra  (w 


where  Rais  the  antenna  radiation  resistance,  given  by  the  real 
part  of  the  antenna  impedance  Z(w). 


Now  from  equations  (15),  it  is  apparent  that  the 
characteristic  behavior  of  the  circuit  is  contained  in  the 
transfer  functions.  For  this  reason,  we  have  calculated  the 
first,  second  and  third  order  transfer  functions.  The  value 
of  the  diode  capacitance  C  was  chosen,  as  by  Kanda  (1980)  to 
be  lO'^^  f.  The  saturation  current  was  also  selected  as  in  Kanda 
to  be  2  X  10'®  amp.  The  diode  exponent  was  calculated  from 

Ct  =  n/qkT 

where  the  diode  ideality  factor  n  was  1,  T  the  temperature  was 
20®  C  (293.15®  K)  and  q  and  k  are  the  electronic  charge  and 
Boltzmann's  constant  respectively.  In  Figure  4,  we  display 
the  result  of  the  calculation  of  the  first  order  transfer 
function  from  equation  (17),  and  the.  scattered  power  from 
equation  (24).  The  power  is  not  the  total  power  at  that 
frequency,  but  rather  what  the  power  would  be  if  all  other 
transfer  functions  were  zero;  in  (24),  we  use  v,(w)  instead  of 
Vp^  The  numerical  data  showed  that  the  first  order  transfer 
function  exhibited  a  maximum  in  modulus  at  525  MHz  and  a 
maximum  of  radiated  power  at  about  550  MHz. 

For  higher  order  transfer  functions,  there  is  a 
representational  problem  that  accompanies  the  display  of 
multidimensional  data.  To  simplify  the  presentation  of  the 
second  order  transfer  function  I^w.,,W2),  we  will  assume  that 
one  of  frequencies  is  fixed  at  the  value  of  27r  x  20  KHz;  this 
audio  frequency  is  chosen  since  it  will  be  used  as  the 
frequency  of  circuit-applied-voltage  source  in  what  follows. 
In  Figure  5,  we  present  the  second  order  transfer  function 
H^w,  ,W2)  with  w;  fixed  and  the  power  radiated  due  to 
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the  second  order  response  only.  The  representation  of  the 
power  is  also  a  bit  confusing:  the  second  order  power  is  also 
a  function  of  2  f requrencies,  one  of  which  is  fixed  at  20  KHz. 
However,  the  power  is  radiated  at  the  frequency  20  KHz  +  w/27r. 
The  discontinuity  in  the  mod  (H)  plot  at  about  100  MHz  is  due 
to  the  transition  from  the  short  dipole  model  to  the 
Schelkunoff  method  of  calcullating  antenna  impedance:  it  is  an 
artifact  and  is  unimportant.  The  maximum  in  radiated  power 
occurs  at  approximately  575  MHz.  The  magnitude  of  the  second 
order  transfer  function  has  2  maxima,  one  at  about  100  MHz, 


and  the  other  at  about  560  MHz. 


Lastly,  we  present  the  third  order  transfer  function 
l^w^  ,W2  ,W3  ).  In  Figure  6,  wi  and  Wj  are  both  fixed  at  20  KHz. 
The  power  in  the  figure  is  radiated  at  frequency  v/2n  +  40 
KHz.  There  are  2  maxima  in  magnitude  of  Hg  ,  at  about  100  MHz 
and  575  MHz  ,  while  the  maximxam  power  is  radiated  at  575  MHz. 
Note  the  similarity  in  the  shape  of  these  curves  and  those  of 
Figure  5. 


A  more  practical  representation  of  the  data  is  by  the 
frequency  terms  in  the  solution  vjC  t )  =v^(  t ) +\^(  t ) +V3(  t )  where 
these  are  given  by  (15).  To  do  this,  we  need  to  specify  the 
source  term  g(t)  =  v(t)/R,.  The  short  circuit  antenna 
current  is  due  to  a  single  plane  wave  incident  broadside  so 
that  at  the  antenna,  the  incident  plane  wave  is 


EoCOS(v(^) 


This  produces  Fourier  coefficients 


i  =  1,2 


(25) 


where  w^  =  wj,  ,  Wj  =  .  Similarly,  the  circuit  voltage  source 

is  also  a  single  sinusoid: 


V  ( t )  =  Vc  cos  ( Ht ) 


which  produces  Fourier  components  of  g(t) 


gj  =  V2Rl  (26) 

where  i  =  3,  4,  and  Wj  =  w^  ,  W4  =  -Wc  . 

Thus,  the  Fourier  series  for  g(t)  has  4  contributions  and  is 
given  by: 


If  we  set  the  circuit  voltage  frequency  w^  to  20  KHz,  an 
audio  frequency,  and  the  incident  plane  wave  frequency  to  500 
MHz,  an  RF  frequency,  the  Volterra  series  solution  will 
produce  frequency  terms  that  are  represented  schematically  in 
Figure  7.  Note  that  some  of  the  terms  are  simply  the  second 
and  third  harmonies.  The  more  interesting  terms  are  the 
mixtures  of  Wg  and  w^  .  Because  Wo»Wc  ,  the  nonlinearity 
results  in  the  appearance  of  sidebands  about  w^  and  2Wq .  The 
pattern  that  emerges  should  be  apparent:  if  we  were  to  carry 

the  Volterra  series  to  fourth  order,  more  sidebands  would 
appear  about  w,,  ,  2Wg  and  they  would  emerge  about  .  The 
usefulness  of  these  sidebands  stems  from  the  antenna  response. 
Although  the  antenna  radiates  negligibly  at  20  KHz,  it 
radiates  strongly  in  the  vicinity  of  500  MHz  and  1  GHz.  Thus, 
the  nonlinearity  results  in  the  incident  RF  plane  wave  being 
scattered  as,  in  effect,  a  carrier  wave  for  the  20  KHz  signal 
in  the  circuit. 

In  Table  1,  we  show  the  results  of  our  simulation  with 
unit  applied  voltage  (i.e.,  1  volt)  at  20  KHz  and  with  a  unit 
incident  plane  wave  (i.e.,  1  volt/m)  at  500  MHz.  The  table 

lists  the  Fourier  components  of  the  solution  at  only  positive 
frequencies;  the  total  solution  for  the  voltage  across  the 
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linear 


f(Hz) 


mod  (Vj  )  arg  (v^  ) 
(volts) 


scattered  power 
(watts) 


20K 

500M 


quadratic 


cubic 


.5 

.  12 


0 

153. 


3.8x10'^* 

7.3x10'® 


40K 

3.9x10'^ 

180. 

3.7x10*3^ 

IG 

0 

500M+20K 

3.0x10* 

7. 8x1  O'* 
2.9xl0* 

51. 

180 

-9.4 

8.7x10*’* 

0 

4.4x10**^ 

.*.“1 

t' 

i 

n 

*  k. 

,  > 

60K 

2.6x10*3 

180. 

8.2x10*3^ 

t 

v\ 

1.5G 

4.6x10'* 

-142. 

5.9x10*** 

» •* 

500M 

20K 

5.8x10''* 

7.7x10*3 

-9.4 

180. 

1.7x10*® 

9. 1x10*3® 

500M+40K 

2.9x10*^ 

-9.4 

4.3x10'*° 

■A 

lGt20K 

5.8x10'* 

51. 

3.4x10-*^ 

1 

Table  1:  Frequency  components  of  simulation  of  Volterra 
series  solution.  Circuit  voltage  source  is  1  volt  at  20  KHz, 
while  the  incident  plane  wave  is  1  volt/m  at  500  MHz. 
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diode  is  constructed  as  the  sum  of  negative  and  positive 
frequency  terms  of  equal  amplitude,  thus  producing  terms  with 
a  time  dependence  of  the  form  cos(wt+0).  Our  focus  here  will 
be  on  the  scattered  or  radiated  power.  The  largest 

scattered  power  is  due  to  the  simple  linear  scattering  of  the 
incident  plane  wave.  Note  that  the  cubic  nonlinearity 
produces  correction  terms  at  both  500  M  and  20  K  that  serve  to 
reduce  the  magnitude  of  the  linear  response  at  these 

frequencies;  this  is  the  compression  term  (Graham  &  Ehrman, 
1973),  a  well  known  linear  effect.  Indeed,  all  nonlinear 
terms  of  odd  order  will  contribute  responses  at  the 
fundamental  frequencies.  For  our  purposes,  the  components  of 
most  interest  are  the  intermodulation  or  mixture  components  at 
500  M  ±  20  K,  1  G  +  20  K  and  500  M  +  40  K.  As  can  be  seen 

from  the  table,  the  largest  radiated  power  occurs  at  500  MHz 

and  the  side  bands  of  +  20  KHz  and  ±  40  KHz.  Since  this 

simulation  is  for  unit  input,  we  can  construct  the  solution 
for  these  frequencies  and  voltage  amplitude  and  plane  wave 
amplitude  Eo  by  using  the  following  simple  proportionality 
relations.  Since  power  depends  on  the  square  of  current,  we 
have  (Pi=  power  for  unit  input); 


f  =  500  M 
f  =  500  M  +  20  K 

f  =  500  M  +  40  K 

f  =  1  G  ±  20  K 

For  example,  if  v^  =  10  volts 

and  Table  1 ; 

f  =  500  M 
f  =  500  M  +  20  K 

f  =  500  M  +  40  K 

f  =  1  G  +  20  K 


P  =  (28) 

P  =  (E^oJ^Pi 
P  =  (E^fPi 
P  =  (I^^P, 

Eo=  1  volt,  we  get  from  (28) 

P  =  7.3  X  id® 

P  =  4.4  X  Id” 

P  =  4.3  X  id® 

P  =  3.4  X  Id” 
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iJote  that  the  component  at  500  M  +  40  K  differs  from  tlie 
linear  term  by  only  a  factor  of  10. 
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We  have  previously  seen  that  the  scattered  power  due  to 
the  second  and  third  order  transfer  functions  is  a  maximum  at 
57‘  'Hz  for  a  circuit  voltage  of  20  KHz.  We  present  the 
result  of  a  simulation  at  these  frequencies  in  Table  2.  A 
comparison  of  the  results  of  Tables  1  and  2  shows  that  the 
intermodulation  components  575  M  +  20  K  and  575  M  +  40  K  are 
larger  than  those  corresponding  to  the  500  M  plane  wave  by 
about  a  factor  of  2.  This  relatively  small  change  is  due  to 
the  relatively  broadband  nature  of  the  antenna. 


f  (Hz  ) 


mod  (Vj,)  arg  (v^  ) 
(volts) 


scattered  power 
(watts ) 


linear 

20K  .5  0 

575M  .125  111. 


3.8x10-22 

3.3x10-5 


quadratic 

40K 

3.9x10“* 

180. 

3.7x10“2^ 

1.  15G 

2.5x10“® 

-35. 

6.7x10*’® 

0 

7.8x10“^ 

180 

0 

575M±20K 

6.0x10*® 

-86. 

7.7x10-’2 

cubic 

60K 

2.6x10“® 

180. 

1.9x10“’^ 

1.725G 

3.6x10“® 

77. 

8.2x10-25 

575M 

1 . 2x1 0“2 

-86. 

3.0x10-® 

20K 

7.7x10“® 

180. 

9.  Ixl0“2® 

575M±40K 

5.9x10“^ 

-86. 

7.5x10-’° 

1. 15G+20K 

4.9x10“^ 

-35. 

2.6xl0“’2 

Table  2;  Frequency  components  of  simulation  of  Volterra 
series  solution.  Circuit  voltage  source  is  1  volt  at  20  KHz 
while  the  incident  plane  wave  is  1  volt/m  at  575  MHz. 


CONCLUSION: 

The  simulations  we  have  performed  have  shown  that  an  RF 
field  is  scattered  by  the  active  diode  circuit  in  such  a 
manner  that  the  RF  frequency  acts  as  a  carrier  for  the  circuit 
signal.  If  the  circuit  has  an  applied  sinusoid  of  20  KHz, 
then,  in  the  absence  of  any  RF  signal,  the  diode  nonlinearity 
results  in  the  voltage  across  the  diode  having  Fourier 
components  of  20  KHz,  40  KHz,  and  higher  harmonics.  With  an 
RF  plane  wave  of  frequency  w^,  incident  on  the  antenna,  the 
resulting  circuit  current  and  scattered  power  has  frequency 
components  of  w^ +  20  KHz,  Wq  ±  40  KHz,  etc.  There  are  also 
components  at  2Wo  +  20  KHz,  2wo +  40  KHz,  etc.  and  the  same  is 
true  for  the  higher  harmonics  of  w^ .  Thus,  each  harmonic  of 
the  RF  signal  acts  as  a  carrier  wave  for  the  circuit  signal, 
in  the  sense  that  the  frequency  composition  of  the  circuit 
signal  is  shifted  from  being  centered  about  zero  to  a  center 
given  by  the  harmonic.  The  relative  magnitudes  of  the  circuit 
signal  are  not  exactly  reproduced  in  the  scattered  waves  at 
the  RF  frequency  and  its  harmonics.  It  should  be  stated  that 
although  we  have  focused  the  discussion  on  the  scattering  by 
the  circuit,  there  are  also  interference  effects  in  the 
circuit  due  to  the  interaction  of  the  circuit  signal  and  the 
RF  signal. 

We  have  performed  calculations  of  a  few  of  these 
scattered  intermodulation  components  using  a  Volterra  series 
approach.  Since  this  is  a  perturbation  method,  the  numerical 
values  derived  can  be  regarded  as  providing  reliable 
information  only  when  our  assumption  of  small  voltage  across 
the  diode  is  appropriate,  or  more  specifically,  from  (8),  when 
tlie  voltage  is  less  than  0.026  volts.  The  data  presented  in 
tlie  tables  violate  this  assumption.  However,  by  appropriate 
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scaling  of  the  unit  responses,  numerical  values  could  be 
easily  generated  which  are  consistent  with  the  perturbation 
approximation.  Nevertheless,  the  formalism  and  numerical 
values  provide  a  semi-quantitative  description  for  the  case 
where  the  approximation  does  not  hold.  A  further  limitation 
of  this  modeling  is  suggested  by  Franceschetti  and  Pinto 
(1980).  These  authors  claim  that  a  model  of  the  diode  which 
ignores  the  nonlinear  junction  capacitive  effects  is  justified 
only  at  frequencies  less  than  1  MHz.  A  more  precise  model 
would  include  these  nonlinear  capacitive  effects  and  would 
carry  out  the  Volterra  series  solution  to  a  higher  order. 
However,  the  computation  of  the  transfer  function  for  orders 
higher  than  five  is  extremely  difficult  and  this  situation 
limits  the  applicability  of  the  Volterra  series  method  to  the 
mildly  nonlinear  case. 
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